Entropic Bounds for the Quantum Marginal Problem 
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The quantum marginal problem asks, given a set of reduced quantum states of a multipartite system, 
whether there exists a joint quantum state consistent with these reduced states. The quantum marginal 
problem is known to be hard to solve in general as it is a variant of the N-representability problem. We 
provide entropic bounds on the number of orthogonal solutions to the quantum marginal problem. 
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The quantum marginal problem (QMP) asks when and 
what joint quantum states of a composite system are com- 
patible with a given set of reduced states. This question has 
its genesis in the marginal problem of classical probability 
theory which is concerned with the existence of a probabil- 
ity density function with given projections onto a set of co- 
ordinate subspaces (see Uj] and references therein). Since 
the foundational work of ||2,|3||, which revealed that almost 
all tripartite quantum states are uniquely determined by 
their two-party reduced states, there has been a great deal 
of progress on variations and generalisations of the QMP 
lliiii&lsKiollilEEOl. At the current time only 
the most fundamental version of the QMP, which asks what 
constraints prevent the existence of a joint state PaiA2---a„ 
given only the single -particle reduced states p^^, is well 
understood. This problem was completely solved by Kly- 
achko in the finite-dimensional setting [SJ, and by Eisert et. 
al. ill in the gaussian setting. 

A general solution to the QMP would have profound and 
revolutionary consequences for physics as it would provide 
a solution to the N-representability problem of quantum 
chemistry ifisifisll . and hence allow us to easily calculate, 
eg., the binding energies of complex molecules. It turns 
out that this is too much to hope for as a general solution 
cannot exist: the A^-representability problem is now un- 
derstood to be too hard to solve, even on a quantum com- 
puter 1IT3I1 . Thus, in order to gain quantitative progress on 
the QMP we must take recourse either to approximate or 
heuristic methods. We take the first approach here: we de- 
rive bounds on the number of orthogonal solutions to the 
general QMP. We illustrate these bounds in the tripartite 
case where we are looking for joint quantum states of a 
tripartite system ABC given two reduced states pab and 
Pbc (even the classical marginal problem is unsolved in 
this case iQl). 

In this Letter we study a generalisation of the QMP, 
namely, how many orthogonal solutions are there to the 
QMP? This problem is of direct relevance not only to quan- 
tum chemistry, but also to condensed matter physics and 
quantum complexity theory because solutions to the QMP 
arise as ground states of locally interacting hamiltonians, 
and the number of orthogonal pure-state solutions is equal 
to the ground-state degeneracy of the hamiltonian. 

We first focus on the case where we are given two re- 



duced states pab and psc; of a tripartite quantum system 
ABC, with local dimensions (Ia, d-B, and dc, and we wish 
to determine how many orthogonal pure states \ipj), j = 
1,2, ... ,171 of ABC are consistent with pab and />bc> 
meaning that tip ( | ip^ ) ( | ) = /O^s and tr ^ ( | ) ( V'j | ) = 
Pbc^ for all j. To approach this problem we suppose that 
m such orthogonal states exist and construct the following 
mixed quantum state of ABC: 



m ^ 



Pabc = 2^—\'^3)Wj\- (1) 

Notice that Pab = trdpABc) and Pbc = trA{pABc)- 

Our next step is to take our putative state Pabc 
and compute its von Neumann entropy, S^bc = 
-iT{pABc^og2{PABc)) = ^og^ifn), which follows be- 
cause the eigenvalues of Pabc are 1/m. 

To complete our derivation we then apply the strong sub- 
additivity inequality for the von Neumann entropy lUol to 
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ABC 



< Sab + S_ 



BC 



ABC, which reads log2(m) 
Sb- Thus we obtain the bound 

^ < + (2) 

Note that Sab, Sbc, and Sb are easy to calculate in terms 
of the initial data as Pb is obtained from pab via partial 
trace. 

By applying the previous argument inductively step we 
can provide a general entropic bound on the number m 
of pure-state solutions to more general instances of the 
QMP: suppose we are given the reduced states pA^Aj+i, 
j = 1, 2, . . . , n — 1 of an n-partite system (the subsys- 
tems Aj need not have the same dimension), then we can 
bound the number m of pure-state solutions to the QMP as 
follows: 

SaiA^-'-A^ ^ SaiA2---A^-i + 5'^„_iA„ ~ Sa^_^ 

1-. „-i (3) 

- + i - 'Y^Sa,, 

where we've repeatedly the strong subadditivity inequality 
to SAiA-i---Ak in each step. Thus, using the same reasoning 
as above, we have that 

71 — 1 n — 1 



•Sa, 



(4) 
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It is clear how to extend this argument to more general set- 
tings. 

The problem of determining the constraints on the ex- 
istence of a joint state p^iAa- a„ given only the single- 
particle reduced states p^. is completely solved. However, 
the number of such constraints is enormous when the di- 
mensions of Aj become large [8]. For example, even for 
the case of two three-level systems there are 387 inequali- 
ties to be checked. For the problem considered here, even 
if these inequalities are satisfied, one still needs to work out 
how many orthogonal pure-state solutions there are. Thus 
it is still desirable to develop bounds on the number of so- 
lutions to the QMP for this simpler case. We can do this 
using the subadditivity of the von Neumann entropy as fol- 
lows. Consider the case where we are given pA and Pb'- 
then, as before, the number m of orthogonal solutions to 
the QMP is bounded by 

log2(m) = Sab <Sa + Sb (5) 

so that m < 2^* 2^^ , and, following the inductive argu- 
ment presented above, we have that 

n 

m < JJ 2^^^ . (6) 

Up to this point we have focussed on the special case 
where we are only looking for pure-state solutions to the 
QMP. However, it could be that there are no pure-state so- 
lutions to the QMP, yet the QMP is still solvable with a 
mixed state. To deal with this we extend our argument to 
the mixed-state case by using another idea, namely, purifi- 
cation [16]: any mixed state p of a quantum system A can 
be realised as the reduced state of a pure state | p) aa' of a 
"doubled" system AA' where A' is a copy of A'. To see 

this write p = X]j=iftl^i)(^il for the spectral decompo- 
sition of p. Then 

d 

\P) = VPiWj) Ai\Uj)A') (7) 

is a purification. If we have purifications 

\p) = Ej=i VWl^j)A(kj)A') and |cr) 

\/^\'^j)A{\vj)*A') of two mixed states p and a 

then 

n 

(pW) = Yl VpJ^i'^j\'^k){{uj\vk)y 

j,k=i (8) 
= tr(Va^) = G{p,a), 

where G{p, a) is a quantity related to the fidelity 
F{p, a) = ii{yrjpa\/p) between p and a. 

We apply this observation in the following way. Sup- 
pose that pabc is a mixed state solution to the QMP where 
we are given pab and psc- Then let \pabca'B'C') be 
the purification Eq. (|7]l of Pabc onto the doubled system 
ABCA'B'C. Suppose that there are m such mixed-state 



solutions with corresponding purifications and sup- 
pose that these purifications are orthogonal. (This orthogo- 
nality implies that the supports of Pabc and Pabc j 7^ ^' 
are orthogonal so that the corresponding solutions Pabc ^o 
the QMP have zero pairwise fidelity: -^(p^^c'' Pabc) ~ 
^jk ) We then apply our main argument to to find 

log2(m) = S ABCA'B'C < 2Sabc 

= 2Sab + 2Sbc ~ 2Sb, 

where we've applied the subadditivity of the von Neumann 
entropy to across the bipartition ABC : A'B'C and 
exploited the fact that Sabc = Sa'B'c- So we learn that 

The extension to the general QMP is clear. 

Note that the arguments presented here don't depend on 
the dimension of the composite system and hence gener- 
alise, in the appropriate limits, to infinite-dimensional sys- 
tems. 

Physically, the bound Eq. ^ says that, in order for there 
to be m solutions to the tripartite QMP the entropy, and 
hence, our ignorance, of the reduced states of AB and BC 
needs to be large enough to suppress the entropy of the in- 
terface system B\ there needs to be enough room to move 
at the interface subsystem to marry up the two reductions 
into a larger consistent state. This is intuitively reasonable 
and provides a physical interpretation for the following 
simple entanglement monogamy result 1, 17]. Suppose we 

are given /9ab = \'^^)ab{'^~\ andpsc = |*~)bc(^"L 
where |*") = ^(|01) - |10)) is the spin-i singlet, then 

there are no pure joint states oabc consistent with these 
reduced density operators iUtII : the bound Eq. (O reads 
m < |, so that m = 0. Our bound Eq. ([Tol l actually 
shows much more as we learn that this situation is robust 
against perturbations; there must be an ball of states around 
\^~)ab and \^^)bc where there are no joint states con- 
sistent with both Pab and Pbc- 

The bounds Eq. (|2l and Eq. ([Tol l are likely to be tight in 
the low-entropy regime, where the entropy of the interface 
system B is larger than that of AB or BC. In this case the 
bounds provide an easy way to prohibit the existence of 
a solution to the QMP. The bound Eq. (|2l) is also likely to 
perform well in the high entropy regime as it reproduces the 
exact result in the completely mixed case Pab = ^/dAdB 
and Pbc = l/dBdc, namely, m = 2'^^+'^^+'^^ . 

In this Letter we have developed an entropic upper bound 
on the number of orthogonal solutions to the quantum 
marginal problem. This bound also provides nontrivial 
constraints on the existence of solutions to the general 
QMP. Using our bound it possible that new bounds for the 
ground-state energy and degeneracy (and hence, the ther- 
modynamic pressure) of local hamiltonians may be devel- 
oped. Additionally, it is likely that our bound will be useful 
in quantum complexity theory where it should be able to 
provide nontrivial bounds on quantum counting problems. 
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